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The truncation or compression of the spectrum of Schmidt values is inherent to the matrix prod¬ 
uct state (MPS) approximation of one-dimensional quantum ground states. We provide a renor¬ 
malization group picture by interpreting this compression as an application of Wilson’s numerical 
renormalization group along the imaginary time direction appearing in the path integral represen¬ 
tation of the state. The location of the physical index is considered as an impurity in the transfer 
matrix and static MPS correlation functions are reinterpreted as dynamical impurity correlations. 
Coarse-graining the transfer matrix is performed using a hybrid variational ansatz based on ma¬ 
trix product operators, combining ideas of MPS and the multi-scale entanglement renormalization 
ansatz. Through numerical comparison with conventional MPS algorithms, we explicitly verify the 
impurity interpretation of MPS compression, as put forward by V. Zauner et al. [New J. Phys. 

IT, 053002 (2015)] for the transverse-field Ising model. Additionally, we motivate the conceptual 
usefulness of endowing MPS with an internal layered structure by studying restricted variational 
subspaces to describe elementary excitations on top of the ground state, which serves to elucidate 
a transparent renormalization group structure ingrained in MPS descriptions of ground states. 

PACS numbers: 03.65.-w, 03.67.Mn, 71.27.-l-a, 75.10.Jm 


I. INTRODUCTION 

In recent years, tensor network states have emerged 
as powerful theoretical and computational tools to in¬ 
vestigate strongly correlated quantum many-body sys¬ 
tems. By focusing on states rather than Hamiltonians, 
tensor networks are purposely designed in order to cap¬ 
ture the entanglement structure inherent in physically 
relevant quantum states. For gapped quantum systems, 
matrix product states (MPS)iI2 are known to faithfully 
represent ground states in one spatial dimensiorP, and 
a plethora of numerical algorithms exist to variationally 
optimize MPS over the manifold of low-energy states of 
local Hamiltonian^^l^. While local quantities can be ap¬ 
proximated to very high precision, long-range behavior 
is not necessarily captured as accurately due to the ex¬ 
ponential decay of correlations in MPS with finite bond 
dimension. This property is of particular importance 
in the context of quantum phase transitions, where it 
has fostered studies of hnite entanglement scaling at 
criticalitjff^®^. 

The multi-scale entanglement renormalization ansatz 
(MERA)P3 is an altogether different kind of tensor net¬ 
work tailored to the description of scale invariant sys¬ 
tems. By introducing unitary disentangling operators, 
its layered structure is able to accommodate a proper 
and sustainable renormalization group (RG) flow along 
the intrinsic scale dimension of the network, even at crit¬ 
icality. Recently, it has been shown that MERA can 
be reinterpreted as stemming from a novel tensor net¬ 
work renormalization (TNR) scheme for coarse-graining 
two-dimensional tensor network^i^lh^. Unlike MPS, the 
MERA incorporates an explicit scale dimension into its 
network structure, which renders its real-space coarse- 
graining properties particularly explicit. Motivated by 


the initial discovery of remarkable similarities between 
holography and MERA^i^, a lot of ongoing worlP^HMl has 
been concerned with clarifying the extent to which it is 
able to realize a true lattice version of the holographic 
principle. 

It is well known that the entanglement content of a 
translation invariant MPS is entirely contained in the 
dominant eigenvector of its transfer matrix. In a recent 
publicatiorpZlj it was observed that the other eigenval¬ 
ues of the MPS transfer matrix contain additional use¬ 
ful information on the elementary excitations and corre¬ 
sponding dispersion relation of the system, providing an 
intriguing connection between the MPS transfer matrix 
and the spectral properties of the Hamiltonian. Put dif¬ 
ferently, there appears to be a highly non-trivial relation¬ 
ship between the excitation spectrum of a local transla¬ 
tion invariant Hamiltonian and the local information and 
static correlations present in its ground state. 

In order to gain a better understanding of the origins 
of the MPS transfer matrix, it was proposed in Ref. [57| 
to treat the physical spin of a MPS as an impurity in 
the two-dimensional tensor network arising from the Eu¬ 
clidean path integral representation of the ground state. 
As depicted in Fig. the transfer matrix of an infi¬ 
nite bond dimension MPS can then be identified with 
the exact quantum transfer matrix at zero temperature. 
By swapping the interpretation of (Euclidean) time and 
space, static correlation functions in the MPS are seen to 
correspond to temporal impurity correlations. The trun¬ 
cated finite-dimensional MPS transfer matrix obtained 
from numerical simulations will therefore contain only 
the degrees of freedom that are relevant to the impu¬ 
rity dynamics up to some infrared cutoff. This per¬ 
spective suggests that the state compression involved 
in variational MPS techniques can be interpreted as an 
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Figure 1. (Color online) Two-dimensional tensor network rep¬ 
resentation of the Euclidean path integral corresponding to 
the ground state of a local one-dimensional translation invari¬ 
ant lattice Hamiltonian H. Horizontal slices are translation 
invariant matrix product operators built from a local tensor 
O (blue circles), and correspond to imaginary time evolution 
with . The MPS ground state |To) is obtained by suc¬ 
cessively applying onto an initial MPS \(p}. Vertical 

columns can then be interpreted as the MPS transfer matrix 
Ta = ® ^*1 where and A' respectively correspond 

to the translation invariant ket and bra MPS matrices. 


where and denote boundary vectors irrelevant 
for bulk properties. The state is hence completely de¬ 
termined by specifying a single tensor A'^ G 
where D and d, respectively, denote the bond dimension 
of the virtual level and the physical dimension of the lo¬ 
cal Hilbert space associated to each lattice site. To each 
MPS we can associate a MPS transfer matrix 

= ^ (2) 

which is a key object in numerical simulations and is used 
to calculate static correlation functions with respect to a 
uniform MPS ground state. We henceforth assume that 
the MPS is normalized such that the transfer matrix has 
a unique largest eigenvalue equal to 1. 

As is well known, a different perspective on the MPS 
transfer matrix can be provided in terms of a two- 
dimensional tensor network associated to the uniform 
MPS representation of the ground state of a local one¬ 
dimensional translation invariant HamiltoniaiPsHUl. For 
clarity, an overview of this construction is given in Fig.j^ 
where we consider the imaginary time evolution 


application of Wilsons numerical renormalization group 
(NRG)P^ along the virtual (imaginary time) dimension 
of the system. 

In this paper, we propose a variational ansatz based on 
matrix product operators (MPOs) to explicitly coarse- 
grain transfer matrices and set out to numerically verify 
the impurity picture proposed in Ref. [27l beyond the an¬ 
alytical results obtained for the XY model in Ref. [HI 
Having a layered decomposition of a MPS ground state 
at our disposal, we benchmark our method on the Ising 
model, propose an ansatz for the structure of MPS fixed 
point reduced density matrices, and study the effect of 
restricting a variational MPS ansatz for elementary ex¬ 
citations to a subspace of variational parameters. We 
furthermore translate our ansatz to the setting of free 
fermions, which allows us to further corroborate our find¬ 
ings by exploiting the additional free fermionic structure. 


II. MPS TRANSFER MATRICES 

Consider a one-dimensional lattice model and a local 
translation invariant Hamiltonian H = hn,n+i re¬ 

stricted to nearest-neighbour interactions. The trans¬ 
lation invariant ground state of this system can be de¬ 
scribed by a uniform MPS 

i'i^[^]) = Evn (1) 
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where \(j)) is an initial state assumed to have non-zero 
overlap with the ground state |^'o)- Due to the lo¬ 
cality of i?, we can split /3 into small imaginary time 
steps (5 and use a Trotter-Suzuki decomposition « 

On to arrive at a translation invariant MPO 

representation of e“^^ (see, e.g.. Ref. jM]). This approxi¬ 
mation introduces a controllable Trotter error depending 
on the order of the decomposition, which renders the 
MPO representation quasi-exact. By grouping tensor 
contractions along imaginary time, we can identify the 
exact MPS transfer matrix with a single column of the 
network. Interpreting the tensor network in Fig. [l]as an 
Euclidean path integral, the exact MPS transfer matrix 
thus coincides with the quantum transfer matrix derived 
from the partition function Z^o = lim^_,.oo Tr (e“^^) at 
zero temperature. Note however that both these ex¬ 
act transfer matrices have exponentially diverging bond 
dimensions, and thus differ from the truncated finite¬ 
dimensional MPS transfer matrix ® 

fined before, which arises in actual numerical simulations. 
Similarly, the exact MPS ground state tensor A® corre¬ 
sponds to a semi-infinite MPO in this picture, and rep¬ 
resents the ground state |4'o) up to some Trotter error, 
of which the truncated finite-dimensional MPS A® is a 
compressed version. Stated in these terms, our goal is to 
shed light on the relationship between the exact MPS A® 
and its compressed version A®. 
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III. COARSE-GRAINING TRANSFER 
MATRICES 

As it is our intention to coarse-grain a translation in¬ 
variant MPO, defined by a local tensor O having on¬ 
site operator dimension d and MPO bond dimension D, 
a natural way to proceed is to devise a coarse-graining 
strategy using MPOs in order to retain the matrix prod¬ 
uct structure. To this end, we introduce an isometric 
coarse-graining MPO denoted by G, which is character¬ 
ized locally by a single five-index tensor 


( 4 ) 


grained tensor 


The tensor g is restricted to be an isometry defined 
by the map g :I(g)I(g)V —^OoV, where I, V, and O, 
respectively, refer to the vector spaces of the incoming 
indices, the virtual indices and the outgoing indices. As 
such, g satisfies g^g — loigiv and gg^ = P, where P is a 
projector onto some subspace of IOIO V. Let us denote 
d = dim(I), X = dim(V) and d' = dim(O). The isometric 
constraints on the local tensor g forces the MPO G as a 
whole to be isometric if > d', or nnitary if d^ = d'. 
Note that there are two different ways to group operator 
and virtual indices of g, leading to two possible “gauge” 
choices, denoted pictorially by the following equations: 



or 



( 5 ) 


Let us now consider a finite set {Gs} of the kind of 
MPO isometries defined above, and label them with a 
discrete scale dimension s = 1,2,, Smax, where each 
Gs is allowed to be different. By acting sequentially with 
each of the {Gs} together with their conjugates {G|} 
on the translation invariant MPO to be coarse-grained, 
we arrive at the MERA-inspired tensor network (see Ap¬ 
pendix]^: 



( 6 ) 


One way to optimize all of these coarse-graining MPOs, 
is to perform a single sweep in scale from bottom to top. 
At each layer s, two tensors Og-i of the previous layer are 
blocked by acting with gs and gl to construct a coarse- 



( 7 ) 


which redistributes entanglement among virtual and op¬ 
erator dimensions. Tracing over the outgoing indices, as 
depicted in Eq. ([^, we can interpret this object as a 
generalized transfer matrix with x^Ps-dimensional fixed 
points ((TL[gs]| and where we have emphasized 

the dependence on the coarse-graining isometry in the 
notation. By locally maximizin^^ the cost functioiP^ 


Gigs 





( 8 ) 


for the isometry gg using a conjugate gradient algorithm 
adapted to unitary manifolda^Zl, we then iteratively up¬ 
date the left and right fixed points of Eq. Q until con¬ 
vergence of Eq. ([^ is attained up to some tolerance. The 
coarse-grained AIPO tensor Og is then obtained by trun¬ 
cating its bond dimension to a fixed number Dg, or up 
to some tolerance, using conventional MPS methodP^^ 
such that 


(9) 


where we have introduced the rank-reducing tensor Wg, 
with denoting the left inverse of Wg such that 

W~^Wg = 1 and WgW~^ projects onto the truncated 
subspace. Note that the new on-site operator dimension 
of the coarse-grained MPO is determined by the output 
dimension d'g of the gate gg. This constitutes the opti¬ 
mization of one layer, and by repeating the above coarse- 
graining procedure until the top level is reached, we ob¬ 
tain a set of effective MPO tensors {Og}, isometric gates 
{gs} and truncation tensors {Wg} for s = 1,2,..., Smax- 
Both the tolerances for the optimization of the cost func¬ 
tion and for the truncations provide sensible measures for 
the errors introduced along the way. 

Although our main motivation in this article is to 
coarse-grain the MPOs appearing naturally along the 
imaginary time direction in the Euclidean path integral 
picture in Eig. there is nothing preventing us from 
applying the above ansatz along the spatial direction 
to renormalize, for instance, a Hamiltonian operator or 
more general transfer matrices arising in two-dimensional 
lattice models. In particular, local Hamiltonian terms 
can be seen to renormalize to a sequence of semi-infinite 
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MPO strings to the left or to the right depending on 
the gauge choice in Eq. ([^, which follows naturally from 
the one-sided causal cone structure that arises due to the 
isometric nature of the coarse-graining gates. 


IV. FROM EUCLIDEAN PATH INTEGRAL TO 
UNIFORM MPS 

The pioneering work of Wilson on the numerical renor¬ 
malization group (NRG) showed that the relevant low- 
energy subspace for an impurity problem could be ob¬ 
tained by applying real-space RG transformation^. Re¬ 
cently, the theory of minimal updates in MERA has re¬ 
lated the success of NRG to the inherent causal cone 
structures arising in MERA^^. The causal cone of a re¬ 
gion had been originally introduced as the part of the 
MERA network that is geometrically connected to and 
able to exert influence on the properties of the state in 
that regiorfi^. It has since been interpreted as the collec¬ 
tion of tensors that needs to be changed in order to ac¬ 
count for a local change of the Hamiltonian in that region, 
which is understood to be sufficient to capture the evolu¬ 
tion of that region under successive coarse-graining trans¬ 
formations while maintaining localitji^. This property 
can be ultimately traced back to the existence of distinct 
energy scales in the Hamiltonian. According to the prin¬ 
ciple of minimal updates, an impurity initially localized 
in space thus remains localized under coarse-graining, 
which leads to a very efficient MERA descript ion o f sys¬ 
tems with boundaries, impurities, or interface^^. The 
relevant degrees of freedom for an impurity are found 
to be exactly those living at the boundary of the causal 
cone. 

Although the MPS ansatz we are about to describe 
bears a resemblance to the impurity MERA construction, 
the underlying motivation is entirely different. We stress 
that the main motivation for our proposal is to be able to 
capture the physics relevant for the physical spin degree 
of freedom of a quantum state, which, as we will argue, 
precisely amounts to extracting the degrees of freedom 
relevant to the spin treated as an impurity in the exact 
quantum transfer matrix. 


A. Construction of MPS ansatz 

Let us now apply our coarse-graining ansatz Eq. 
to the vertical transfer matrix Ta obtained from the 
imaginary time evolution network in Fig. Note that 
the physical dimension and virtual dimension in Fig. 
now correspond, respectively, to the bond dimension and 
the operator dimension of the coarse-graining MPO. We 
furthermore implicitly assume that Ta = i.e. the 
transfer matrix is understood to originate from some Eu¬ 
clidean rotated effective local Hamiltonian m In gen¬ 
eral, there is no reason to expect the Hamiltonian H to 


be related to the physical Hamiltonian H involved in the 
MPO description of the imaginary time evolution . 

To arrive at a uniform MPS, we will treat the physi¬ 
cal spin connecting the exact MPS representations of ket 
and bra at r = 0 as an impurity in Ta- The arbitrary 
location of the physical index in imaginary time does not 
a priori introduce an inhomogeneity in Ta- However, for 
expectation values of local operators different from the 
identity, the privileged nature of the physical spin be¬ 
comes manifest, and an operator insertion at r = 0 leads 
to a modified transfer matrix Tj’ = used 

in calculating MPS expectation values. As such, a static 
correlation function between two operators separated by 
n sites in the physical system corresponds to a “temporal 
correlator” of the impurity between two operators sepa¬ 
rated by n steps of evolution with Ta = where H 
again denotes not the physical Hamiltonian H but the 
Euclidean rotated effective local Hamiltonian. We will 
thus construct a truncated MPS representation A® by 
extracting the relevant degrees of freedom for the “Eu¬ 
clidean dynamics” of the impurity from the exact quan¬ 
tum transfer matrix Ta- 

In what follows, we will focus on the upper semi-infinite 
MPO describing the ket part |4'[A]) of the transfer ma¬ 
trix Ta, as depicted in Fig.[^ and refer to the location of 
the physical spin as the “impurity” regardless of whether 
there is actually an operator inserted at r = 0. To explic¬ 
itly construct the approximated ket tensor A® associated 
to the uniform MPS |'k[A]) in Fig. we first apply our 
ansatz to the infinite quantum transfer matrix, and op¬ 
timize for Smax layers as if no impurity were present. We 
then insert the impurity at r = 0 and draw its causal 
cone C in Fig. [^a), where the inside of the causal cone 
stretches out to the left and to the right, and the degrees 
of freedom affecting the impurity at r = 0 arise from con¬ 
tracting the network outside the causal cone. By ignoring 
the bottom half of the network, we can identify the part 
of the coarse-graining network which acts on the semi¬ 
infinite ket part of the MPO. In retrospect, the impurity 
naturally suggests which isometric restriction in Eq. ([^ 
to impose on the coarse-graining tensors during the op¬ 
timization, as the entire tensor network surrounding the 
original semi-infinite transfer matrix can be interpreted 
as a projector onto the dominant eigenvector subspace 
relevant to the impurity, see Fig. [^b). 

It is clear that only one index of the MPO coarse- 
graining tensors crosses the boundary of the causal cone 
for each layer. Together, these legs constitute the sites of 
an effective lattice system defined along the bound¬ 
ary of the causal cone, which is called the Wilson chain. 
Sites along this chain are labeled by the layer index 
s = 1, 2,..., Smax, where site s contains an effective 
renormalized description, for d' = d, of 2® sites of the 
original lattice located roughly at a distance 2® away from 
the impurity. Moving along the Wilson chain thus cor¬ 
responds to changing scale and moving away from the 
impurity. Next, we insert approximate resolutions of the 
identity at the lowest layer in Fig. [^c) using the rank- 




5 


scale (s) 



Figure 2. (Color online) Uniform MPS from coarse-graining the exact transfer matrix. Location of the physical spin, considered 
as an impurity at r = 0, is denoted with a red cross, (a) Bonds crossing the one-sided causal cone C of the impurity at r = 0 
identify degrees of freedom relevant to the physical spin, which live on the Wilson chain defined along the boundary of 
C. The tensor network outside C renormalizes the exact transfer matrix Ta into a new transfer matrix Tw along . (b) 
Contraction of the tensor network outside C renormalizes the ket, projecting the exact MPS A* onto the subspace defined by the 
relevant degrees of freedom along . (c) Insertion of approximate resolutions of the identity 1 ~ WiW^^. (d) Coarse-graining 
the semi-infinite MPO string to the next layer leaves a single Wi tensor behind. (e)(f) Repeating this procedure leads to the 
construction of the Wilson MPO. (g) After introducing an infrared cutoff by capping the Wilson MPO with the fixed point of 
the uppermost coarse-grained tensor, we retrieve a finite-dimensional MPS A*. 


reducing tensors WiW^^ obtained during the optimiza¬ 
tion to reduce the bond dimension in the imaginary time 
direction. As shown in Fig. |^d), the semi-infinite MPO 
string is pushed to the next coarse-graining layer, leav¬ 
ing a single tensor Wi behind. Repeating this procedure 
for the next layer in Fig. [^e)-(f) leads to the emergence 
of an inhomogeneous MPO along the imaginary time di¬ 
rection, which we will refer to as the Wilson MPO, and 
which effectively amounts to a sequence of coarse-grained 
Trotter step^. 


By introducing the fixed point of the uppermost 
coarse-grained tensor as an infrared cutoff in Fig.j^g), we 
arrive at a finite-dimensional approximate Wilson-based 
MPS with an internal layered structure resulting in the 
uniform MPS tensor A*. The reason for this particu¬ 
lar cutoff strategy can be traced back to interpreting the 
transfer matrix Ta as a thermal state (or more generally 
a mixed state) with exponentially decaying correlations. 
We then expect the coarse graining network to be able to 
disentangle this state into a product state using a finite 
number of layers, similar as in the case of ground states 
of gapped Hamiltonians, such that, after a finite number 


of layers, we can pictorially denote the flow by 



If we want the coarse-grained MPO tensor at the top 
to have bond dimension = I, the optimal choice 

is given by the eigenvectors corresponding to the largest 
eigenvalue of Eq. Q. For a critical MPO, we expect 
scale invariance to direct the coarse graining process to 
a fixed point that needs to be iterated forever. In that 
case, introducing the product state fixed point as an in¬ 
frared cutoff explicitly breaks the critical properties and 
can as such only be considered an approximation, whose 
quality can be assessed by evaluating the energy with 
respect to the resulting state. Note that all variational, 
finite bond dimension MPS approximations necessarily 
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contain an implicit infrared cutoff, and that different cut¬ 
off implementations might yield different MPS tensors 
approximating the same gapless state, a point to which 
we return in Sec. IV Al 

In the following two sections, we consider two imme¬ 
diate applications made possible by the structure of our 
ansatz. 


B. Structure of MPS fixed point reduced density 
matrices 


From the layered structure of our MPS decomposition, 
it is straightforward to associate a layered decomposition 
to the zero-dimensional fixed point reduced density ma¬ 
trices pl and pa of the MPS transfer matrix Eq. ([2) as 
well. To see this, let us apply a truncation procedure 
to the transfer matrix 7^ constructed from our ansatz, 
instead of working solely on the level of the ket |'I'[A]). 
By sequentially grouping indices across bra and ket for 
the truncation step at each layer, we eventually arrive at 
the norm 



where the grey triangles denote the Wilson MPOs for 
ket and bra, and the red triangles are now invertible 
rank-reducing tensors (left) and Xg (right), for 

s = 1,2,..., Sniax, that can be obtained from lifting the 
impurity from bottom to top and truncating the effec¬ 
tive transfer matrix at each layer s by considering it as a 
MPS, 




Note how we can approximately determine the fixed 
points of the transfer matrix in a sequential manner with¬ 
out explicitly calculating them, and how only small bond 
dimensions are involved in the intermediate truncation 
processes. A potential application of this insight lies 
in the problem of efficiently contracting two-dimensional 


tensor networks, e.g. overlaps of projected entangled pair 
state^^, where the one of the main bottlenecks in vari¬ 
ational ground state optimizations can be traced back 
to calculating high-dimensional MPS fixed points of the 
MPO transfer matrices needed for the calculation of envi¬ 
ronment tensors. The feedback of information from high- 
to low-energy in our layered decomposition of these hxed 
points might lead to a substantial reduction in compu¬ 
tational resources, thus allowing the ^proximation of 
otherwise intractable MPS fixed points^. 

Another application concerns scale invariant theories, 
for which we expect to find a recursive relation, directly 
from the layered MPS decomposition, whose eigendecom- 
position ought to contain the scaling operators and scal¬ 
ing dimensions of the critical theory that we are approx¬ 
imating. Suppose we have identified a fixed point trun¬ 
cation tensor W* and associated X* in the layered de¬ 
composition of our MPS ansatz, characterizing the scale 
invariant fixed-point behavior. We can then propose a 
“radial” superoperator 



(13) 


whose eigenvalues and eigenvectors are expected to be 
related to the scaling dimensions and operators of the 
underlying conformal field theory. Forgetting about the 
Wilson-based MPS origin of the ansatz, we can also 
build concentric tensor networks from Eq. (13), which en¬ 
force scale invariance and can b e regarded as variational 
ansatze in their own rightP^^. Similar superoperators 
for ext ractin g scaling fields were of course conceived in 
MERAI^^I^ and TNRP^, but not in a MPS setting. A 
partial approach to extract scaling information from a 
MPS was investigated in Refs. ITOHI^ 


C. Restricted variational subspaces for excitations 

Given the internal layered structure of the uniform 
MPS in Fig. [^g), it is tempting to ask how excitations 
can be interpreted in this framework. Following Wilson’s 
RG interpretation of the Kondo impurity probleirP^, we 
expect low-energy excitations to live near the top of the 
Wilson MPO if the layered MPS decomposition obtained 
from our ansatz is to be interpreted as a true renormal¬ 
ization group scale. Note that our ansatz retains trans¬ 
lation invariance in space as we effectively construct a 
uniform MPS, and that the scale dimension in our case 
refers to coarse-grained imaginary time. We want to 
emphasize that, even though there is no explicit spatial 
coarse-graining taking place (as the number of sites re¬ 
mains constant), spatial correlations do arise naturally 
in our picture as we coarse-grain imaginary time to grow 
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Figure 3. (Color online) Example of a restricted variational 
MPS excitation ansatz (_Bs)) at intermediate layer index 
s, constructed to allow for targeting the variational parame¬ 
ters X residing in an interval of layers [s, Smax] (green trian¬ 
gles), while the tensors corresponding to layers [1, s) are held 
fixed. 


the virtual bond dimension of the MPS. Despite the ab¬ 
sence of an explicit coarse-graining in space, the succes¬ 
sive MPO layers labeled by s are non-unitary and can 
therefore shrink the effective Hilbert space at every step. 
The idea is now to expose this effective reduction of the 
Hilbert space at higher renormalization scales by per¬ 
turbing the tensors within the layered structure of the 
MPS and observing to which part of the spectrum of the 
Hamiltonian this particular perturbation gives access. 

In particular, exploiting the translation invariance 
which we still have at our disposal, we can use the varia¬ 
tional MPS ansatz for localized excitations developed in 
Refs. SSI and SHI given by 

|<l>(^)(i?)) = 

— \m<n / \m'>n / 

(14) 

Comparing with Eq. ([^, it is clear that the excitation 
ansatz is constructed on top of the ground state wave- 
function, as we change a single ground state tensor 
into a tensor H® and take a momentum superposition 
of this localized perturbation. All variational freedom is 
contained within the tensor H®, and the variational opti¬ 
mization of the Rayleigh-Ritz quotient 


{4^'>{b)\h\4^\b)) 
s (<i>(^)(i?)|4^)(H)) 


(15) 


gives rise to a generalize eigenvalue problem. In order for 
this generalized eigenvalue problem to be well-defined, it 
will be necessary to project out so-called null-modes, i.e. 
(almost) zero eigenvalues of the effective normalization 
matrix in the right hand side of the generalized eigenvalue 
equation. 

Let us now recast Eq. (14) pictorially and substitute 
our layered MPS decomposition of the MPS ground state 
to arrive at the variational ansatz depicted in Fig. i By 
fixing part of H® to be equal to A® using the layered de¬ 
composition, we can now design a restricted excitation 


ansatz for s = 1, 2,..., Smax, which only al¬ 

lows for variations with respect to an interval [s, Smax] 
of tensors in the Wilson MPO description of the full ex¬ 
citation ansatz. The restricted tensor R®(x) = Ss(a:^) is 
defined in terms of a linear mapping 

Xa'X«' £DxdxD 

which glues the variational parameters x to the surround¬ 
ing MPS ground state tensors A® G ^DxdxD _ p^j. 
ferent layers s, this tensor parametrizes different sub¬ 
spaces of the full variational spac^^ spanned by the 
states |4>p"^^(i3)) with unrestricted H®. Starting from 
^Smax’ '''®hore only the variational degrees of freedom in 
the top tensor are considered, we can gradually take all 
layers [1,..., Smax] into account until the ansatz culmi¬ 
nates into the full = B\. 

Note that the role of our restricted ansatz is not so 
much in improving the efficiency of numerically calcu¬ 
lating excitations within the framework established in 
Ref.HHl but in providing a novel and conceptually intrigu¬ 
ing interpretation of excitations within MPS in a way 
that explicitly tries to resolve the different energy scales 
present in the MPS ground state tensor. In Sec. |VB[ we 
will verify that the number of null-modes as function of 
layer index s and momentum p is intimately tied to the 
effective reduction of the Hilbert space at higher renor¬ 
malization scales. 


V. NUMERICAL RESULTS 

To assess the validity of our ansatz, we have performed 
numerical simulations using the MPS ansatz introduced 
in Sec. II V Al and discuss the transfer matrix and Schmidt 
spectra of the resulting MPS ground states. We calcu¬ 
late dispersion relations using the setup of Sec. |IV C| to 
illustrate that variational subspaces corresponding to the 
upper layers of our MPS ansatz are effectively restricted 
to the low-energy part of the spectrum. Additionally, we 
translate our method for compressing transfer matrices 
to the setting of free fermions, where the additionally 
imposed structure on the numerics allows us to highlight 
and further corroborate our findings. 


A. MPS ansatz from Wilson MPO 

Consider the quantum Ising model in a transverse field, 
which can be defined in terms of the Hamiltonian 

i7 = -E<^f+i-AE'"^ (17) 

i i 

where A determines the strength of the applied magnetic 
field. One can easily show how to construct a trans¬ 
lation invariant MPO representatiorP^ of the Trotter- 
Suzuki decomposition of , to which we apply our 
ansatz along the imaginary time direction. Comparing 
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Figure 4. (Color online) Comparison of transfer matrix (TM) 
spectra (left) and Schmidt spectra (right) of the transverse 
Ising model Eq. ( |17[ ) between our MPS ansatz 14*[A]) based on 
the Wilson MPO (blue circles) and a variational MPS ground 
state |4'[^]) (red crosses) for A = 1.1 (gapped) and D = 16. 
All MPO simulations were performed using a Trotter step S = 
0.001 and local dimensions ds,Xs,d's € {2,4} of the coarse- 
graining tensors with Smax = 14 layers. 


the ground state energy density of the MPS resulting 
from our ansatz with the exact solution, we find rel¬ 
ative energy errors AE = {E — Eq)/Eq ~ 10“® both 
for the gapped (A ^ 1) and gapless (A = 1) phases us¬ 
ing very modest computational resources. Having ob¬ 
tained a uniform MPS from the Wilson MPO, we can 
further quantify the accuracy of our ansatz by study¬ 
ing the spectrum of the transfer matrix Eq. ([^, and the 
Schmidt values of the MPS, which can be retrieved from 
the spectrum of the dominant eigenvector of the transfer 
matrix Eq. ([2|. In Fig. i] and Fig. we show a compar¬ 
ison of low-lying transfer matrix spectra obtained with 
respectively our coarse-graining method and a variation- 
ally optimized MPS using the time-dependent variational 
principle. 

For the gapped case depicted in Fig. an excellent 
match is found for both the TM eigenvalues A„ and the 
Schmidt coefficients A^ as soon as the number of coarse- 
graining layers is chosen sufficiently large, which demon¬ 
strates the behavior of the infrared cutoff as explained in 
Sec. ITV^ 


For the gapless case in Fig. the infrared cutoff can 
be recognized as a very particular initial state in the ef¬ 
fective imaginary time evolution induced by the Wilson 
MPO, potentially spoiling the relevant infrared data for 
critical models. From the impurity point of view, there 
is a strong interaction between the different sites in the 
Wilson chain corresponding to degrees of freedom 
living at different energy scales. The infrared cutoff in¬ 
troduced after a finite number of layers at one end of the 
Wilson chain, has a strong feedback on the physics of the 
impurity, which is living at the other end of the Wilson 
chain. As argued in Sec. |IV A different implementations 
of the implicit infrared cutoff for gapless states, yield 
different variational MPS approximations, which has re¬ 
cently been discussed in Ref. UHl 

To illustrate this explicitly, we also studied the effect 
of further truncating the Wilson based MPS to some spe- 




Figure 5. (Color online) Comparison of transfer matrix (TM) 
spectra and Schmidt spectra of the transverse Ising model 
Eq. ( |17[ ) between our MPS ansatz |4'[A]) based on the Wil¬ 
son MPO (black circles) and a variational MPS ground state 
in'[A]) (crosses) with bond dimension D = 12 (red), D = 24 
(blue) and D = 50 (black) for A = 1 (gapless). All MPO sim¬ 
ulations were performed using a Trotter step 5 = 0.001 and 
local dimensions ds,Xa,d's G {2, 4} of the coarse-graining ten¬ 
sors with Smax = 17 layers, and we also investigated the effect 
of further truncating the MPO based ansatz to specific bond 
dimensions D = 12 (red circles) and D = 24 (blue circles) 
based on the Schmidt coefficients. 


cific bond dimensions using the standard MPS recipe 
(j.e. throwing away the smallest Schmidt coefficients). 
While this has no distinguishable effect on the remain¬ 
ing Schmidt coefficients of Fig. ib) , the effect on the 
nondominant part of the spectrum of the transfer matrix 
itself is significant, as can be observed in Fig. [^a). Nev¬ 
ertheless, when comparing to a variational MPS with a 
particular bond dimension, we obtain a qualitative agree¬ 
ment that matches the first few dominant Schmidt coeffi¬ 
cients (see also Fig. for analogous entanglement spec¬ 
tra obtained for critical free fermions). 


B. MPS excitation ansatz from Wilson MPO 

Using the restricted variational subspace ansatz de¬ 
fined in Fig. we have calculated variational approxi¬ 
mations to the dispersion relation of the transverse Ising 
model Eq. 0 for all layers intervals [s, Smax], for s = 
1,2,..., Smax, of the Wilson MPO. As depicted in Fig. 
sweeping across layer intervals gradually allows less layers 
to contribute to the variational optimization. We observe 
that restricting the variational degrees of freedom to the 
top layers limits the momentum range of the ansatz to 
such an extent that, within the accuracy of the ground 
state approximation itself, momentum states correspond¬ 
ing to high-energy excitations can no longer be captured. 
Diluting exponentially on their way down the network, 
momentum states corresponding to high-energy excita¬ 
tions are seen to yield states with a norm that effectively 
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layer interval [s, Smax = 14] 


0 

momentum p 


Figure 6. (Color online) Variational approximations to the 
quantum Ising dispersion relation e{p) = 2-y/l + — 2Acosp 

(solid black line) in function of the layer index s for A = 
1-1, Smax = 14 layers and bond dimension D = 16 using 
the restricted variational MPS excitation ansatz of Fig. 
Momenta near the minimum p = 0 of the dispersion relation 
are seen to be fully captured by restricting the variational 
ansatz to the top tensor(s) of the Wilson MPO, in contrast to 
momenta corresponding to high-energy excitations, which are 
observed to correspond to null states with a norm effectively 
approaching zero. By varying the range of the variational 
degrees of freedom in scale space from infrared towards the 
ultraviolet, we recover the high-energy excitations. 


approaches zero, leading to ill-defined generalized eigen¬ 
value problems Eq. (15) as a function of layer index s 
and momentum p. This observation of null-modes sug¬ 
gests an effective reduction of the Hilbert space of the 
momentum states corresponding to high-energy excita¬ 
tions at coarser renormalization group scales. 


C. Free fermion ansatz 

We now develop a very similar approach to coarse- 
graining transfer matrices using free fermions. For the 
Ising model in Eq. we can exploit the fact that it 

is solvable by mapping it onto a system of free fermions 
to gain independent evidence in support of our ansatJ^. 
The additional free fermionic structure will furthermore 
allow us to study explicitly what happens to the transfer 
matrix on every coarse-graining level at the level of the 
ferminionic modes, which is impossible to extract from 
the MPO construction. To that end, we emplo y the 
modified version of the free fermionic MERAE21S2I to ef¬ 
fectively construct the network of coarse-graining MPOs 
appearing in Eq. Q. For technical details regarding the 
construction, as well as further details, we refer to Ap¬ 
pendix]^ 

The renormalized transfer matrix T® is depicted in 
Fig. I^a). At each renormalization scale s it can be di¬ 


agonalized as 


T® = exp 


T c]L Cn 


(18) 


where Cm are free fermionic annihilation operators, and 
the spectrum is fully determined by the single parti¬ 
cle “energies” e^. The coarse-graining MPOs are con¬ 
structed layer by layer so as to properly describe the 
dominant low-energy part of the spectrum of the virtual 
Hamiltonian H. As such, the free fermionic nature of 
the problem allows us to keep track of the spectrum of 
the transfer matrix at each coarse-graining step, as pre¬ 
sented in Fig. [^a,b) for non-critical and critical systems 
respectively. 

We then use the coarse-grained MPOs obtained above 
to sequentially compress the transfer matrix, as depicted 
in Fig. m- This allows us to observe how the spec¬ 
trum of the compressed transfer matrix 77®jp is gradually 
emerging with growing s, where s = 0 marks the original 
transfer matrix and s = Smax is a fully compressed one, 
see Fig. I^a) for comparison. We plot the single parti¬ 
cle energies equivalent to Eq. (181, in Fig. |^c,d) 

both for non-critical and critical systems. The spectra 
consist of a continuous part which can be attributed to 
still-to-be-renormalized low-energy part of the spectrum, 
and discrete high-energy modes corresponding to the im¬ 
purity degrees of freedom. There are two such modes 
emerging at each new layer (for = 2^, one comes from 
the ket and the other one from the bra). We can therefore 
argue that each fermionic mode, supported on the Wil¬ 
son chain, represents a different, exponentially shrinking 
momentum-window (recall that the original problem was 


s s 




Figure 7. (Color online) (a) Coarse-grained transfer matrix 
and (b) partially compressed transfer matrix, obtained by ap¬ 
plication of s layers of coarse-graining MPO. The same isom¬ 
etry is used here for the top part (i.e. ket) and the bottom 
part (i.e. bra) of the original MPS. 
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Figure 8. (Color online) Single particle spectrum of (a,b) the 
coarse-grained transfer matrix T" in Fig. [7|a), and of (c,d) 
the impurity transfer matrix Ti^p in Fig.TTTb), at different 
layers s. In (c,d) the approximately continuous part of the 
spectrum is drawn with a solid line, while pluses mark the 
discrete high-energy modes which are localized at, and can be 
associated with, the impurity degrees of freedom. ResultJ^ 
for non-critical A = 1.1 (a,c) and critical A = 1 (b,d). The 
initial transfer matrix was obtained by applying 4096 layers 
of MPOs with bond dimensions equivalent to Xs = 2^ and 

ds=d's= 2*. 


translationally invariant), analogous to Wilson’s renor¬ 
malization group picture of the impurity problem. 

The spectrum of the compressed transfer matrix is then 
effectively given by 

where e{k) is the dispersion relation of the original trans¬ 
fer matrix Ta , and fc™P ~ A™ is a logarithmic discretiza¬ 
tion of momenta resulting from the RG scheme, at least 
up to some infrared cutoff related to the finite length- 
scale inevitably appearing in the problem. 

Most importantly, exactly the same structure of the 
transfer matrix spectrum was observed resulting from the 
conventional MPS truncation procedur^^. Note, how¬ 
ever, that here the parameter A (A = V2 in Fig. is 
a constant which depends on the geometry of the tensor 
network through the bond dimensions ds,Xs, d[, while for 
the standard truncatiorP^ it depends non-trivially both 
on the bond dimension of the MPS and correlation length 
appearing in the problem (possibly as the result of finite 
entanglement scaling for the critical system). As such, 
we can expect that the geometry of the tensor networks 
limits the precision of the ansatz in a similar way to us¬ 
ing a finite bond dimension in MPS, which we discuss in 
more detail in Appendix [Aj 

Finally, in order to provide direct evidence that our 
ansatz allows for effective compression of the transfer 
matrix, we show that it is reproducing the entanglement 
spectrum with good accuracy. To that end, the reduced 
density matrix of the dominant eigenvector of can 



index m index m 


Figure 9. (Color online) Dominant part of the entanglement 
spectrum for increasing scale s. The spectrum is preserved 
with good accuracy. Results for (a) A = 1.1 and (b) A = 1^. 
Bond dimensions equivalent to Xs = 2^ and (a) ds = d'g = 2^, 
(b) ds = d's = 2®. The initial transfer matrix was obtained 
by applying 4096 layers of MPOs. See the text for further 
discussion. 


be expressed aP^ 


P = ^ 


-2j2SUUn 


( 20 ) 


where fm are fermionic annihilation operators, Z 
is the normalization factor, and the entanglement 
spectrum is determined by 6^. Consequently, the 
Schmidt coefficients, up to normalization, are given by 
.... In Fig. we plot the domi¬ 
nant part of the spectrum for growing s, both (a) for the 
non-critical system A = 1.1 and (b) the critical one A = 1. 
Notice that for the non-critical case s = Smax = 10 cor¬ 
responds to truncating down to 10 fermionic modes, and 
the plot shows 7 dominant modes in the entanglement 
spectrum, proving that the compression is quite effec¬ 
tive. For the critical case we stop at s = 8, just before 
applying the top tensor. As encountered in Sec. |V A[ 
simple top tensors limit the number of modes describing 
the state too quickly, which can significantly affect the re¬ 
sulting spectrum. We can further increase the precision 
by changing the geometry of the tensor network, either 
by increasing x, or equivalently, by compressing slower, 
which amounts to increasing d' for fixed d, for which we 
refer to Appendix [A} 


VI. CONCLUSIONS 

In this paper, we have shown how the state compres¬ 
sion inherent to variational MPS techniques can be inter¬ 
preted as resulting from a renormalization group proce¬ 
dure applied to the Euclidean path integral description 
of a quantum system. By treating the location of the 
physical spin as an impurity, we were able to construct a 
uniform MPS representation that takes into account the 
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degrees of freedom relevant to the impurity. The MPO 
structure of the coarse-graining ansatz led to a natural 
proposal for the structure of the MPS fixed point re¬ 
duced density matrices. Furthermore, we explicitly re¬ 
lated the different layers in the decomposition with en¬ 
ergy scales by studying a restricted variational ansatz for 
excitations, which shows that perturbations at the high 
layers only give access to elementary excitations with a 
momentum near the minimum of the dispersion relation, 
whereas perturbations with other momenta get washed 
out by the subsequent layers below it. We also formu¬ 
lated a free fermion version of our coarse-graining ansatz 
that is in agreement with the results of the MPO ansatz. 

Having arrived at a conceptually suggestive picture of 
the renormalization group structure inherent to matrix 
product states, we can look at the possibility of con¬ 
ceptual advantages in calculating correlation functions, 
scattering matrices, and other quantities. Further study 
of the behavior of excitations at this boundary between 
MPS and entanglement renormalization might yield in¬ 
sight into how to develop a proper excitation ansatz for 
MERA, which is even more relevant in light of its pur¬ 
ported relation to holography. Finally, as our Wilson- 
based MPS furthermore suggests that the burden of en¬ 
tanglement can be shifted to manageable correlations be¬ 
tween energy scales, it would be interesting to explore the 
possibility of continuum generalizations of our ansatz in 
terms of the continuous MERA^S^, and to study its appli¬ 
cability to the numerical optimization of two-dimensional 
quantum latti ce sy stems using projected entangled pair 
states (PEPSj^S^. 
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Appendix A: Free fermion construction 


1. Transfer matrix 


Suzuki-Trotter decomposition of an operator e 


V = 


Vi = exp 




(Al) 

(A2) 

(A3) 


for some smalP^^ time-step S. Here, V represents a single 
row in the two-dimensional network in Fig. It appears 
naturally as the transfer matrix operator in the solution 
of the two-dimensional classical Ising model and as such 
was extensively studied, see e.g. the Review lATl Relevant 
for us, V has a simple representation in terms of an MPO 
with bond-dimension 2 (see e.g. Ref. IMjl . which allows us 
to directly obtain the full transfer matrix single 

column in Fig.[^, of the same form as the V above. 


ta = wy^w2wy^, 

Wi = exp ^ATi ^ 

W 2 = exp ^K 2 ^ T-f : 


(A4) 

(A5) 

(A6) 


where are standard Pauli matrices acting on the 
virtual degrees of freedom labeled with I = 1,2,...2L, 
where L is the number of times V was applied onto the 
initial state (the ground state is obtained in the limit of 
L —>• 00 ). The physical degree of freedom is localized at 
the bond between sites L and L -|- 1. Finally, the pa¬ 
rameters Ki ^2 can be found as ATi = — | logtanh((5) and 
K 2 = — \ logtanh((5A). 

The transfer matrix Ta can be diagonalized in a stan¬ 
dard wajE3 by mapping onto a system of free fermions 
with the Jordan-Wigner transformation = ia 2 n-io- 2 n, 
Tn = a2n-i rinKn Convenience, we intro¬ 

duce Major ana fermions am n = 1,2,3, ...4L, which 
are Hermitian by construction and satisfy the canoni¬ 
cal anti-commutation relations {an, am} = 2Sn,m- The 
transfer matrix can then be diagonalized as 7a = 

exp Y^m=i ^mii> 2 m-ib 2 m^ , where the index m can be 
identified with momentum and is the dispersion rela¬ 
tion. The Majorana fermions bm = O^^^n^n are described 
by the orthogonal matrix . Since the transfer matrix 
has effectively open boundary conditions, some care is 
needed during diagonalization. To that end, we follow 
the procedure outlined in Ref. [Ml obtaining and 
numerically for some large, fixed value of L. 


2. Coarse-graining procedure 

In order to construct the tensor network representation 

of the ground state of the quantum Ising model in Eq. Our main goal is the construction of (the equivalent of) 

0. as depicted in Fig. we start with the second-order the coarse-graining MPO in Eq. Q. To that end, the 
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transfer matrix, which is a Hermitian, positively defined 
and Gaussian operator, can be described uniquely (up to 
a normalization) using the correlation matrix 


(^Ta _ 

*^771,n 


TtTa 


(A7) 


We use the free fermionic MERA ansatz, adapting the 
construction described iiP to our problem. The coarse- 
graining MPO is decomposed here as 



P 


U 




]Q 


(A8) 


where e; are single particle energies of the correspond¬ 
ing transfer matrix. Optimizing the MERA then boils 
down to bringing the correlation matrix as 

close as possible to the canonical form, which allows for 
the identification of the high-energy modes that need to 
be truncated during the coarse-graining. We note that 
the truncation procedure is equivalent to tracing out the 
high-energy modes in the transfer matrix, which could, 
in principle, introduce some bias when compared to the 
(more correct) projection on their ground state. We ob¬ 
serve however that the procedure is working well, see 
Fig.|8]:a,b). 


3. Compressing the transfer matrix 


consisting of one disentangler U and one isometry W in 
the language of MERA. The gates U and W, which are 
in principle position and layer dependent (we skip the 
indices to simplify notation), are parametrized in terms 
of SO(2P) matrices describing the local canonical trans¬ 
formation of Majorana fermions. They are shifted with 
respect to each other by one site (two Majorana modes), 
fixing the equivalent of the bond dimension in the coarse- 
graining MPO to X = 2^, with suitable boundary con¬ 
ditions for the tensors at the ends of the chain. The 
isometries ^ where 

2Q 2P-2Q 

y- = diag(lX'^^Xl,0;or^^TM), (A9) 

and y+ = 1 — Y~ are diagonal 2P x 2P matrices selecting 
the first 2Q and the last 2P — 2Q Majorana fermions, 
respectively. A single layer of the coarse-graining MPOs 
in Eq. Q is in this picture a direct sum of disentanglers 
and isometries in that layer, and we will mark it as W^. 
The coarse-grained correlation matrix is obtained from 
the previous one as 

c® = (AlO) 

with representing the coarse-grained trans¬ 

fer matrix depicted in Fig. [^a). 

We optimize the orthogonal matrices U and W in each 
layer, starting with the bottom layer and progressing 
to the top, using the standard optimization strategy for 
MERA^^ The cost function is given by 

f,{U,W) = Tr (W+l'C'^-^W+W^), (All) 

which allows for a local optimization. In the above. 


As indicated in the main text, we use the MERA gen¬ 
erated above to construct the isometry used to compress 
the transfer matrix. We argue that the impurity (physical 
spin) can be well described using the degrees of freedom 
living on the boundary of its light cone, see Fig. To 
that end, we construct the isometry C/j^p which partially 
compresses the transfer matrix up to the scale s. 


(A14) 


where full compression is achieved for s = s^ax- The 
partially compressed 7(^p, described by the correlation 
matrix Cfj^jp, is obtained by acting with on the cor¬ 
relation matrix as indicated in Fig. j^b). The same 
t/j^p is used for ket and bra of the original MPS matrix, 
so only about half of the tensors generated with MERA 
are used here in the end. 

By bringing the correlation matrix Cfj^^p into its canon¬ 
ical form, we can find the spectrum of the corresponding 
transfer matrix as 



Tlmn = exp 


-E' 


imp,s t 


Cm Cr) 


(A15) 


in its diagonal base of Dirac fermions Cm- The spec¬ 
trum is plotted in Fig. [^c,d) and discussed further in 
the main text. The entanglement spectrum, Eq. (201, 


which is given by the reduced density matrix of the dom¬ 
inant eigenvector ( 
well, see e.g. Ref. 


inant eigenvector of 77^p, can be calculated from Cf^^^p as 


4. Comparison of impurity picture with standard 
truncation 


Notice that the canonical form of the correlation matrix, 
following a suitable canonical transformation, is 



1 

itanhe;/2 


—ztanh ei/2\ 

1 r 


(A13) 


Finally, we can now directly compare the isometry 
17()jjp obtained with MERA, with the standard trunca¬ 
tion procedure described by the isometry C/s„d, a 2L x 2% 
matrix where 2^ is the bond dimension of the truncated 
MPS. The matrix Usvd is obtained by calculating the 
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Figure 10. Dominant part of the entanglement spectrum and the overlap between the isometries Usvd and (see the text) 
for the first couple of coarse-graining layers s. Results for (a) A = 1.1 and (b) A = 1^. The initial transfer matrix was obtained 
by applying 4096 layers of MPOs. Comparison of different rates of compression in one layer, set by the ratio of QjP, see Eq. 
(A8L This corresponds to the bond dimensions Xs = “2^, ds = 2^^^ and d' = 2^ of the coarse-graining MPO. 


reduced density operator of the MPS on a half-infinite 
chain (our initial A’' in Fig. Ilj) and finding dominant 
modes in its diagonal basis^Sl. To that end, we look at 
the singular values of which directly show how 

well the dominant modes in the entanglement spectrum 
are preserved during the compression. 


The results for the gapped and critical cases are plot¬ 
ted in Figs. 10 (a) and (b) respectively, together with the 
resulting entanglement spectra. Notice that each mode in 
the entanglement spectrum corresponds to two Majorana 
modes in term of the isometries U. We show the results 
for the first couple of la yers, up to s = 5, and various 
ratios of Q/P, see Eq. (A8). This ratio sets an effec¬ 
tive light cone, and determines what fraction of the large 
energy modes of the original transfer matrix is renormal¬ 
ized during each coarse-graining step. Notice that from 


that perspective, one step for Q/P = 2/8 is equivalent 
to two steps for Q/P = 4/8 and almost five steps for 
Q/P = 6/8. At the same time, the number of modes 
describing the impurity resulting from each step is the 
same. The trade-off between the rate of compression and 
the accuracy can be readily seen. As already described in 
the main text, we expect that this is directly equivalent 
to obtaining MPS with given bond dimension D using 
the standard procedure. For instance, for Q/P = 2/8 
only the first few modes are described accurately corre¬ 
sponding to smaller D, while for Q/P = 6/8 all the domi¬ 
nant (as plotted) modes are reproduced corresponding to 
large D (we should be able to obtain the same result by 
increasing y). The choice oi Q/P = 4/8 (i.e., d = d') 
seems to be a good compromise between the two. 
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Appendix B: Isometric MPO versus MERA 

In the main text, we have proposed to characterize an 
isometric matrix product operator (MPO) locally by in¬ 
troducing a five-legged tensor g, which can be interpreted 
as a map g : —>■ OoV, where I, V, and O, respec¬ 

tively, refer to the vector spaces of the incoming indices, 
the virtual indices and the outgoing indices. As a vari¬ 
ational set of states, the ansatz Eq. § , constructed of 
layers of isometric MPOs, includes the MERA. Explicitly, 
the ansatz Eq. ([^ encompasses MERA in the sense that 
when the local isometry g is, for example, interpreted as 
a disentangler u and an isometry w of a binary MERA, 



the internal bond connecting the disentangler and the 
isometry is in no way restricted hy d = dim(I). Note, 
however, that there is no local gauge transformation act¬ 
ing purely on the virtual level of the coarse-graining 
MPO, which can transform between Eqs. (B11 and (B21. 
Intuitively, we expect that the virtual dimensions of the 
MPO structure within every layer result in a coarse grain¬ 
ing scheme that is quasi-local rather than strictly local 
as in the MERA). This difference is due to the different 
causal cone structures. The causal region of the MPO 
ansatz is potentially larger than that of MERA as it, 
at most, extends towards infinity in one of both direc¬ 
tions. The MERA causal cone, which is strictly local, 
is included in this extended causal cone. The more gen¬ 
eral MPO ansatz still enables us to renormalize trans¬ 
lation invariant MPOs into translation invariant MPOs 
and renders the identification of the Wilson chain very 
straightforward. 
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